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Abstract 

Ice  accreted  on  structures  from  treezing  rain  causes  both  increased  vertical 
loads  and  increased  wind  loads,  due  to  the  larger  projected  area  of  the  structure. 
Structural  failures  initiafed  by  ice  loads  frequently  cause  millions  of  dollars  of 
damage  fo  overhead  power  and  communication  lines,  towers,  and  other  ice- 
sensitive  structures.  There  is  little  information  on  ice  loads  to  use  in  the  design 
of  these  structures,  so  freezing-rain  models  have  been  developed  for  use  with 
weather  measurements  to  determine  the  severity  of  accreted  ice  loads  from 
historical  data.  This  report  describes  a  detailed  heat-balance  ice  accretion 
model,  including  the  important  heat  fluxes  in  freezing  rain  and  allowing  the 
accretion  of  runoffwater  in  the  form  of  icicles.  Italso  presents  a  simple  algorithm 
tor  calculating  the  ice  load  on  components  with  different  diameters  and  cross 
sections.  Collision  efficiency  in  freezing  rain  and  the  calculation  of  the  wind-on- 
ice  load  are  also  discussed.  Model  results  are  compared  with  the  ice  load 
measured  during  a  recent  freezing  rain  storm,  and  to  each  other  using  45  years 
of  weather  data  from  Des  Moines,  Iowa. 


Coven  In  this  1 994  freezIng-raln  storm  In  the  Southeast 
Ice  accretions  on  power  lines  and  trees  caused 
over$l  00  million  In  damage.  (Photo  courtesy  of 
Entergy.) 


Forconversion  of  SI  units  to  non-SI  units  of  measurement  consult  ASTM  Standard 
E380-93a,  Standard  Practice  for  Use  of  the  International  System  of  Units, 
published  by  the  American  Society  for  Testing  and  Materials,  1 916  Race  St., 
Philadelphia,  Pa.  19103. 
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NOMENCLATURE 

Cp  specific  heat  of  air  at  constant  pressure  = 
1.006  J/g°C 

Cpc  heat  capacity  of  the  wire  (J/g  °C) 

Cw  specific  heat  of  water  =  4.22  J/g  °C 
Cd  drag  coefficient 

d  di'op  diameter  (pm) 

D  diameter 

Dq  wire-ice  diameter  (cm) 

Di  icicle  diameter  (cm) 

dd  diameter  of  protoicicle  =  0.5  cm 

dewT  dew-point  temperature  (°C) 
dt  time  interval  (s) 

eo  saturation  vapor  pressure  of  air  over  water 

at  T  =  0°C  =  6.1121  (1.0007  +  3.46  x  10"^  Pa) 
mbar 

ej  saturation  vapor  pressure  of  air  over  water 

at  temperature  T  =  6.1121  (1.0007  +  3.46  x 
10"^  Pa)  ell7.502T/(240.97  +  T)]  mbar 
/  fraction  of  impinging  precipitation  that 
freezes  to  the  wire 

fij  fraction  of  runoff  water  that  forms  a 
protoicicle 

fs  fraction  of  rimoff  water  that  freezes  to  sides 

of  an  existing  icicle 

g  acceleration  of  gravity  =  980  cm^/  s 
Gr^  Grashof  number  (based  on  characteristic 
length  L)  =  g(0  -  T)L^/(273.15  +  T)Va^ 
h  heat  transfer  coefficient  (W / °C) 

hjxi  vapor  transfer  coefficient  (m/s) 

I  electric  current  (amp) 

fca  thermal  conductivity  of  air  = 

4.186  X  10”^  (573  +  1.8  T)  W/m  “C 
L  length  scale  for  nondimensional  numbers 

Lbar  average  icicle  length  =  Li/ 45  cm 
Le  latent  heat  of  vaporization  for  water  =  2501 

j/g 

Li  total  icicle  length  (cm) 

Lf  latent  heat  of  fusion  for  water  =  334  J/g 

M  ice  mass  on  wire/length  (g/m) 

Mi  icicle  mass / length  of  wire  (g/m) 
nta  molecular  weight  of  dry  air  (g/ mol) 
me  mass  of  evaporated  water /length  of  wire 
(g/m) 

mw  molecular  weight  of  water  vapor  =  18.016 
g/mol 

Ni/l  Nusselt  number  (based  on  characteristic 
length  L)  =  /zLAa 

P  precipitation  rate  (mm /hr) 

Pa  atmospheric  pressure  (mbar) 

Pr  Prandtl  number  =  |Xa^^p/fca 

Ps  solar  power/length  (W /m) 


Qc  convective  heat  flux  (W  / m^) 

Qe  evaporative  heat  flux  (W /m^) 

Qf  latent  heat  flux  (W/m  ) 

Ql  longwave  heat  flux  (W / m^) 

Qk  heat  flux  from  droplet  kinetic  energy 
(W/m^) 

Qr  resistive  heating  flux  (W/m^) 

Qs  shortwave  heat  flux  (W /m^) 

Qv  viscous  heat  flux  (W / m^) 

Qw  heat  flux  from  warming  impinging  precipi¬ 
tation  (W/m^) 

R  gas  constant  =  8.31441  J/mol  K 
Rai^  Rayleigh  number  (based  on  characteristic 
length  L)  =  GrLPr 

Rc  resistance/length  of  conductor  (ohm/ m) 

Rei^  Reynolds  number  (based  on  characteristic 
length  L)  =  LV/Va 

Req  equivalent  uniform  radial  ice  thickness  (cm) 
Rh  relative  humidity  =  edewl/^T 
Ty  overall  recovery  factor  for  viscous  heating 
=  0.79 

S  perimeter  of  cross  section 

Sc  Schmidt  number  =  Va/K^ 

S/zL  Sherwood  number  (based  on  characteristic 

length  L)  = 

Sr  solar  radiation  flux  (W/m^) 

T  air  temperature  (°C) 

V  wind  speed  (m/s) 

l^load  wind  load /length  = 

Pay^(Dc  +  Di  Li/100)/2  N/m 

V  drop  velocity  (m/s) 

W  liquid  water  content  (g/m^) 
w  flux  of  water  (g/ m^  s) 

a  angle  from  horizontal 

At  time  step  =  1  hour 

8  ratio  of  the  molecular  weights  of  water 
vapor  and  air  =  nta 

Km  molecular  diffusivity  of  water  vapor  in  air 
=  0.211  [(T  +  273.15)/273.15]^-^^  (1013/Pfl) 
cm^/s 

Pa  dynamic  viscosity  of  air  = 

1.71  X 10“^  +  5.2  X  10~^T  g/cm  s 
Va  kinematic  viscosity  of  air  =  Pa^Pa  cm^/ s 
Pa  density  of  air  =  0.348  x  10“^  Pa/{T  +  273.15) 
g/cm 

Pc  wire  density  (g/cm'^) 

Pi  assumed  density  of  glaze  ice  =  0.9  g/cm^ 
Po  density  of  water  =  1.0  g/cm^ 

PwT  saturated  vapor  density  at  temperature  T 
(g/cm^) 

a  Stefan-Boltzmann  constant  =  5.67  x  10  ° 
W/m^ 
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Ice  Accretion  in  Freezing  Rain 


KATHLEEN  F.  JONES 


1.  INTRODUCTION 

Freezing  rain  occurs  throughout  the  United 
States.  It  may  cause  hazardous  conditions  for  pe¬ 
destrians  and  cars  as  the  rain  freezes  to  sidewalks 
and  streets.  The  rain  may  also  freeze  to  structures. 
Open  structures  that  are  not  designed  to  with¬ 
stand  the  accreted  ice  load  and  the  increased  wind 
load,  due  to  the  greater  projected  area  of  the  struc¬ 
ture,  may  fail.  Recent  severe  freezing  rain  storms 
in  the  Midwest  and  the  South  that  caused  struc¬ 
tural  failures  initiated  by  ice  and  wind  are  de¬ 
scribed  in  Storm  Data  (NOAA  1959-1995:  Mar. 
1990;  Oct,  Nov.  1991;  Feb.,  Mar.  1994).  Both  trans- 
mission-line  design  guidelines  (lEC  1990,  ASCE 
1991)  and  tower  standards  (EIA/TIA 1991)  recom¬ 
mend  designing  for  accreted  ice,  but  provide  little 
information  on  the  ice  load  that  should  be  used  in 
design. 

There  have  been  a  few  nationwide  and  regional 
studies  to  establish  the  severity  of  freezing  rain  in 
this  country.  The  only  systematic  measurements 
of  accreted  ice  thicknesses  in  the  United  States  are 
described  by  Bennett  (1959),  who  presents  data 
from  nine  years  of  measurements  in  the  1920s  and 
1930s.  Tattelman  and  Gringorten  (1973)  used 
Bennett's  data  and  information  from  Storm  Data 
and  its  predecessors  to  determine  ice  thicknesses 
from  freezing  rain  storms  in  the  United  States  over 
a  50-year  period.  They  then  calculated  extreme  ice 
loads  for  return  periods  of  up  to  100  years  in  each 
of  seven  regions  of  the  country.  Vilcans  and 
Burnham  (1989)  mapped  the  number  of  hours 
with  freezing  rain  in  the  24  years  from  1965  to  1988 
at  110  airport  weather  stations  in  the  eastern  half  of 
the  U.S.  using  archived  weather  data.  MRI  (1977) 
estimated  ice  loads  on  transmission  lines  for  25-, 
50-,  and  100-year  return  periods  in  Washington, 
Oregon,  Idaho,  and  part  of  Montana  using  weather 


data  with  their  proprietary  freezing  rain,  snow, 
and  rime  ice  accretion  model.  The  Southeast  Re¬ 
gional  Climate  Center  (1993)  determined  the  aver¬ 
age  number  of  freezing-rain  events  per  year  and 
the  average  event  duration  in  Virginia,  North  and 
South  Carolina,  Georgia,  Alabama,  and  Florida 
using  archived  weather  information. 

Hourly  weather  information  from  himdreds  of 
weather  stations  across  the  country  has  been  com¬ 
piled  at  the  National  Climate  Data  Center  (NCDC) 
and  its  military  counterpart,  the  Environmental 
Technical  Applications  Center  (ETAC).  Weather 
data  are  available  in  magnetic  tape  format,  begin¬ 
ning  in  the  1940s  for  some  stations.  These  hourly 
measurements  of  the  weather  conditions,  includ¬ 
ing  precipitation  rate,  air  temperature,  dew-point 
temperature,  wind  speed,  atmospheric  pressure, 
and  solar  radiation,  can  be  used  with  an  ice  accre¬ 
tion  model  to  hindcast  the  amount  of  ice  accreted 
in  past  freezing-rain  storms. 

This  report  describes  two  freezing-rain  models 
I  developed  to  use  the  historical  hourly  weather 
data.  The  first,  described  in  sections  4  and  5,  is  a 
detailed  heat-balance  model  that  determines  the 
amount  of  the  impinging  rain  that  freezes  either 
directly  to  a  structure  or  as  icicles.  The  second  is  a 
simple  flux  model  that  requires  only  precipitation 
and  wind-speed  measurements  to  determine  the 
amoimt  of  accreted  ice,  assuming  that  all  the  im¬ 
pinging  water  freezes.  This  simple  model,  pre¬ 
sented  in  section  6,  gives  physical  insight  into  the 
factors  governing  the  ice  load  on  structural  ele¬ 
ments  with  different  cross  sections.  Both  models 
are  incorporated  in  the  Microsoft  FORTRAN77 
program  ZRAIN.*^  The  weather  parameters  used 
in  the  two  models  are  described  in  section  2.  Ac- 


*Available  by  request  from 
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tual  and  modeled  ice  accretion  shapes  in  freezing 
rain  are  discussed  in  section  3.  The  factors  affect¬ 
ing  the  collision  efficiency  of  freezing  rain  drops, 
in  comparison  with  cloud  droplets,  are  explained 
in  section  7.  The  calculation  of  the  wind  load  on  an 
ice-covered  structure  using  modeled  ice  loads  is 
presented  in  section  8.  Accreted  ice  thicknesses 
and  masses  predicted  by  the  two  models  are  com¬ 
pared  in  section  9. 

2.  HOURLY  WEATHER  DATA 

Weather  elements  required  in  the  detailed  heat- 
balance  model  are 

Present  weather  code 
Wind  speed 
Precipitation  rate 
Air  temperature 
Dew-point  temperature 
Solar  radiation 
Atmospheric  pressure. 

Only  present  weather  code,  precipitation  rate,  and 
wind  speed  are  used  in  the  simple  model.  Conve¬ 
nient,  but  not  consistent,  imits  are  used  in  ZRAIN: 
wind  speed  is  in  m/s,  precipitation  rate  in  mm/hr, 
solar  radiation  in  W/m^,  atmospheric  pressure  in 
mbar,  and  air  temperature  and  dew-point  tempera¬ 
ture  in  °C.  The  dew-point  temperature,  along  with 
the  air  temperature,  is  used  to  determine  the  rela¬ 
tive  humidity  of  the  air: 

Rh  =  edewx/^T ,  (1) 

where  ej;  and  ^dewT  are  the  saturation  vapor  pres¬ 
sures  of  air  over  water  at  temperatures  T  and  dewT, 
respectively.  The  diameter  of  the  structural  ele¬ 
ment  on  which  the  ice  load  is  to  be  determined  is 
specified  in  centimeters.  Appropriate  factors  are 
incorporated  in  the  model  to  convert  to  a  consis¬ 
tent  set  of  imits. 

Neither  model  predicts  freezing  rain.  They  both 
rely  on  the  weather  observer  having  indicated  by 
the  present  weather  code  when  freezing  rain  occur¬ 
red.  A  number  of  present  weather  codes  are  used 
to  record  various  intensities  of  freezing  drizzle, 
freezing  rain,  and  mixed  precipitation.  ZRAIN 
assumes  that  the  user  has  screened  the  weather 
data,  extracted  the  freezing-rain  events,  and  sum¬ 
marized  the  occurrence  of  precipitation  in  each 
hour  as  freezing  rain  or  drizzle  (Z),  mixed  precipi¬ 
tation  (+),  rain  (R),  or  snow  (S)  by  examining  the 


four  present  weather  code  fields  in  the  meteoro¬ 
logical  records.  The  present  weather  code  is  as¬ 
sumed  to  apply  for  the  entire  hour  it  is  reported. 
Only  freezing  rain  or  drizzle,  rain  with  the  air 
temperature  below  freezing,  or  mixed  precipita¬ 
tion  freeze  to  the  structure.  In  particular,  the  mod¬ 
els  do  not  attempt  to  calculate  the  accretion  of 
snow. 

ZRAIN  also  requires  the  history  of  the  anemom¬ 
eter  height  above  ground  at  the  weather  station. 
The  wind  speed  at  a  specified  height  above  ground 
is  obtained  from  the  measured  wind  speed,  using 
the  anemometer  height  (ft)  and  assuming  a  1/7 
power  law  and  exposure  C  (ASCE  1994).  Typi¬ 
cally,  anemometer  heights  in  the  United  States 
have  varied  between  15  and  100  ft  above  ground, 
so  the  correction  of  the  wind  speed  to  10  m  (the 
worldwide  default  anemometer  height)  or  30  m 
(typical  height  of  a  transmission-line  wire)  maybe 
significant  (Fig.  1).  Anemometer  height  histories 
for  first-order  weather  stations  are  available  in 
Local  Climatological  Data  Annual  Summary  (NOAA, 
yearly),  and  for  military  stations  from  ETAC  in 
electronic  format. 

Solar  radiation  data  are  used  to  determine  the 
solar  heat  flux  at  the  accretion  surface.  Typically, 
there  is  little  solar  radiation  during  freezing  rain 
because  of  the  complete  cloud  cover  and  the  rela¬ 
tively  low  solar  angle  in  winter.  During  the  day, 
solar  radiation  fluxes  in  these  conditions  range  up 
to  about  200  W /m^,  with  values  of  about  80  W/m 
more  common.  If  solar  radiation  information  is 
not  available,  then  it  can  be  assumed  that  there  is 
no  incoming  solar  radiation.  This  is  true  at  night, 
but  this  assumption  will  reduce  the  total  heat  flux 
from  the  actual  amount  during  the  day  and  may 
slightly  increase  the  hindcast  ice  load. 

Atmospheric  pressure  is  used  in  calculating  the 
density  pa  and  saturation  vapor  pressure  ej  of  air 
and  the  diffusivity  of  water  vapor  in  air  Km-  The 
dependence  of  these  parameters  on  atmospheric 
pressure  is  not  great.  Fora  50-mbar  pressure  change 
from  1000  mbar,  pa  and  Km  change  by  5%  and  ejhy 
only  0.02%.  If  atmospheric  pressure  data  are  not 
available,  then  the  pressure  appropriate  for  the 
elevation  of  the  site  may  be  assumed,  with  little 
effect  on  the  hindcast  ice  load. 

Most  weather  elements  are  measured  every 
hour  and  are  reported  as  an  instantaneous  value, 
rather  than  as  an  average  over  the  hour.  In  the 
heat-balance  model  I  assume  that  temperature, 
dew-point  temperature,  atmospheric  pressure,  and 
wind  speed  vary  linearly  between  the  reported 
values.  In  each  hour  there  are  six  10-min  time 
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Figure  1.  Wind  profile  for  a  10-m/s  wind  at  10  tn  above  ground 
assuming  the  If?  power  law  (ASCE  1994), 


steps.  In  the  first  time  step  (e.g.,  1400-1410)  these 
weather  elements  take  on  the  values  reported  for 
that  hour  (1400),  and  in  the  last  time  step  (1450- 
1500)  they  are  equal  to  the  values  reported  for  the 
next  hour  (1500).  In  each  intermediate  time  step, 
one-fifth  of  the  difference  between  the  hourly  val¬ 
ues  is  added  to  the  value  for  the  previous  time  step. 
This  scheme  weights  the  reported  values:  each 
applies  for  20  minutes,  and  the  interpolated  values 
apply  for  10  minutes.  Precipitation  rate  and  solar 
radiation  are  reported  as  a  total  for  the  hour,  and 
I  assume  them  to  be  constant  during  each  hour.  In 
the  simple  model,  the  wind  speed  is  also  assumed 
to  be  constant  for  each  hour. 

A  missing  precipitation  rate  for  an  hour  in 
which  precipitation  is  reported  by  the  present 
weather  code  is  assigned  0.01  mm/hr  in  ZRAIN. 
This  gives  continuity  to  the  freezing-rain  event 
without  arbitrarily  adding  significant  mass  to  the 
accretion.  If  wind  speed,  temperature,  dew-point 
temperature,  or  air  pressure  is  missing  in  an  hour 
it  is  assumed  to  be  the  same  as  in  the  previous 
hour. 

The  models  contain  no  procedure  for  melting 
accreted  ice.  I  assume,  in  extracting  freezing-rain 
events  from  the  weather  data,  that  accreted  ice 
melts  as  soon  as,  and  only  when,  the  air  tempera¬ 


ture  rises  above  freezing  after  the  freezing  rain 
ends.  The  occurrence  of  T  >  0°C  after  freezing  rain 
ends  thus  defines  the  end  of  an  event.  Many  freez¬ 
ing-rain  storms  end  with  warming  temperatures; 
however,  after  a  well-documented  ice  storm  in 
Hanover,  N.H.,  the  temperature  remained  cold, 
and  accreted  ice  remained  on  wires  for  more  than 
30  days  (Ackley  and  Itagaki  1970). 

3.  ICE  ACCRETION  SHAPE 

In  the  models,  the  structure  on  which  the  ice  is 
accreting  is  a  horizontal  circular  cylinder,  which  is 
referred  to  as  a  wire  in  the  following.  Extrapola¬ 
tion  of  the  models'  results  to  other  cylinder  cross 
sections  is  discussed  in  section  6. 

Ice  accreted  on  horizontal  cylinders  (wires  and 
tree  branches)  during  freezing-rain  storms  in  New 
York,  Arkansas,  and  Vermont  is  shown  in  Figure 
2  and  on  the  cover.  The  shape  of  the  ice  accretions 
varies  from  a  crescent  on  one  side  of  a  branch  (Fig. 
2d),  to  a  round  accretion  with  a  few  icicles  (Fig.  2a), 
to  a  heavily  icicled  barbed- wire  fence  (Fig.  2b)  and 
power  line  wire  (cover),  to  an  elongated  accretion 
shape,  apparently  fused  icicles,  on  a  wire.  The 
shape  of  accreted  ice  depends  on  the  diameter  of 
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the  wire,  the  rate  of  precipitation,  the  wind  speed, 
and  the  rate  at  which  the  impinging  water  freezes. 
The  models,  however,  do  not  attempt  to  deter¬ 
mine  the  ice  accretion  shape,  but  instead  assume 
simple  shapes  that,  I  hope,  adequately  represent 
reality. 

In  the  simple  flux  model,  the  ice  accretes  with  a 
uniform  radial  thickness  around  the  circumfer¬ 
ence  of  the  wire.  This  simple  shape  is  consistent 
with  the  level  of  detail  in  the  model. 

In  the  heat-balance  model,  the  precipitation 
that  freezes  immediately  accretes  with  a  uniform 
radial  thickness.  The  runoff  water  that  does  not 
freeze  immediately  is  allowed  to  freeze  as  icicles. 
The  "icicle"  in  this  model  is  a  single,  vertically 
oriented,  circular  cylinder  of  ice  per  meter  of  wire 
that  can  be  understood  as  the  actual  icicles  stacked 
up  end  to  end.  As  runoff  water  freezes  to  the  sides 
of  the  icicle,  its  diameter  increases  from  an  initial 
assumed  diameter  of  0.5  cm  (Maeno  et  al.  1994). 
Runoff  water  that  does  not  freeze  to  the  sides  of  the 
icicle  is  assumed  to  flow  down  to  its  tip  where  it 
may  freeze,  in  successive  0.5  cm-diameter  half 
spheres,  and  increase  the  length  of  the  icicle.  Water 
that  does  not  freeze  there  is  assumed  to  drip  off.  At 
the  end  of  each  time  step,  the  total  icicle  mass  is 
distributed  evenly  over  the  new  icicle  length  to 
form  a  uniform  circular  cylinder.  For  each  10-min 
time  step  the  diameter  of  Ihe  wire-ice  accretion  Dc, 
and  the  diameter  Di  and  length  Li  of  the  icicle  are 
determined  as  initial  conditions  for  the  next  time 
step. 


4.  HEAT  BALANCE  IN  FREEZING  RAIN 

The  primary  process  in  determining  the  amount 
of  the  available  precipitation  that  freezes  to  a 
structure  is  the  flux  of  heat  from  the  accretion 
surface.  Heat  must  be  removed  from  the  rain  drops 
that  collide  with  a  structure  for  them  to  freeze.  In 
calculating  the  heat  fluxes  in  the  heat-balance 
model,  the  wire-ice  accretion  is  assumed  to  be  at 
0°C.  This  is  usually  a  good  assumption  in  freezing 
rain  because,  in  the  typical  near-freezing  air  tem¬ 
peratures,  relatively  low  winds,  and  high  water 
flux,  not  all  the  impinging  precipitation  freezes. 

The  flux  of  precipitation  is  computed  using  the 
relationship  from  Best  (1949)  to  determine  the 
liquid  water  content  of  the  air  as  a  function  of  the 
precipitation  rate, 

W  =  0.067P°-^^  (2) 

where  W  is  in  g/m^  and  P  is  in  mm/hr.  Best 


reviews  relationships  between  precipitation  rate 
and  liquid  water  content  that  both  he  and  other 
researchers,  including  Marshall  and  Palmer  (1948), 
had  determined.  I  chose  his  formula  rather  than 
the  Marshall-Palmer  formula  (W  =  0.072P^'^^)  for 
ZRAIN  because  it  is  the  average  of  all  the  W(P) 
relationships,  while  Marshall-Palmer's  formula 
gives  comparatively  high  values  of  W,  In  ZRAIN, 
the  raindrops  move  horizontally  at  the  wind  speed, 
so  the  flux  of  precipitation  w  is  the  vector  sum  of 
the  vertical  and  horizontal  water  fluxes,  Ppo  and 
Wy,  respectively,  converted  to  consistent  units: 


w{g/m^s) 


■1/2 


+  {WVf 


(3) 


Nine  heat-flux  terms  were  considered  for  inclu¬ 
sion  in  the  heat-balance  model  to  determine  the 
fraction  of  the  available  precipitation  that  freezes. 
They  are: 


•  Qw  (cooling  flux  required  to  warm 
the  incoming  water  to  0°C)  = 

-Cy^zv{0  -  T) 

•  Qc  (convective  cooling  flux)  = 

-nh{0  -  T) 

•  Qe  (evaporative  cooling  flux)  = 
-lo\hniLe[eo/273.15  -  Rhej/iT 
+  273.15)]mw/P 

•  Qi  (longwave  radiative  cooling  flux) 

=  -7ia[273.15^  -  (T  +  273.15)^] 

•  Qs  (shortwave  radiation  heat  flux)  = 

nSr/2  (4) 

•  Qf  (flux  of  heat  released  by  freezing 
all  impinging  precipitation)  =  LfW 

•  Qv  (viscous  heating  flux)  = 

10-^nhr^V^/{2cp) 

•  Qk  (heat  flux  from  droplet  kinetic 
energy) = 

lQ-^w{V^  +  [Ppo/(3.6W)]^)/2 

•  Qr  (resistive  heating  flux  from  cur¬ 
rent  flowing  through  a  conductor)  = 

10¥  Rc/Dc. 


The  heat  fluxes  are  in  watts  per  unit  length  of  the 
wire  per  wire-ice  diameter.  Qo  Qe/  Ql/  Qs/  and  Qv 
all  include  a  factor  of  k  because  they  are  for  the 
outside  surface  of  the  wire-ice  accretion,  which 
has  an  area  per  unit  length  of  n  times  the  diameter. 
The  magnitudes  of  these  fluxes  (except  Qr)  are 
compared  in  Figure  3  for  ranges  of  temperatures, 
wind  speeds,  precipitation  rates,  relative  humid- 
ties,  and  diffuse  solar  radiation  typical  for  freez¬ 
ing-rain  events.  The  individual  heat-flux  terms  are 
described  below. 
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e.  Relative  Humidity 


Figures.  Comparison  of  heat  flux  terms  in  typical  freezing-rain  conditions.  Standard  conditions  (black  symbols)  for 
the  parameters  areV  -S  m/s,  T  =  -3°C,F  =  3  mm/hr,  Rh  =  0.90,  Sr  =  100  W/m^,  Dc  =  3  cm.  In  each  plot  one  or  two 
of  these  parameters  is  varied  about  the  standard  conditions  as  follows:  a)  0  <V  <  10  m/s,  b)  -6  <T  <  -0.5°C, 
c)  0.5  <  P  <  10  mm/hr,  d)  0.5  <  Dc  <  20  cm,  e)  0.50  <  Rh  <  2.00  and  0<Sr<  200  W/m^.  Minimum  convective  and 
evaporative  fluxes  in  a)  are  provided  by  free  convection.  Key:  e  =  evaporative  cooling,  c  =  convective  cooling,  I  = 
longwave  cooling,  w  =  coolingfrom  warming  water,  s  -  shortwave  solar  heating, f=  heat  of  fusion,  v = viscous  heating, 
k  =  kinetic  heating. 


Cooling  flux  required  to  warm  impinging 
precipitation 

In  the  heat-balance  model,  the  latent  heat  of 
fusion  given  off  by  the  fraction  of  precipitation 
that  freezes  to  the  wire  warms  all  the  impinging 
precipitation  to  0°C.  Assuming  the  raindrops  are 
at  the  air  temperature  T,  Qw  is  the  cooling  flux 
necessary  to  raise  their  temperature  to  0°C. 

Convective  cooling 

One  of  the  primary  processes  that  cools  the 
wire-ice  accretion  is  the  loss  of  heat  from  the  cold 
air  blowing  past  the  structure.  When  there  is  no 
wind,  heat  is  removed  at  a  slower  rate  by  free 
convection  controlled  by  the  temperature  differ¬ 
ence  between  the  wire-ice  accretion  and  the  air. 
The  Nusselt  number  Nwl  is  used  to  define  the  rate 
of  heat  transfer  in  both  forced  and  free  convection 
in  external  flow: 

Nwl  =  hl/k^,  (5) 

where  h  is  the  heat  transfer  coefficient,  L  is  the 
appropriate  length  scale,  and  k^  is  the  thermal 
conductivity  of  the  air.  In  forced  convection,  Nu'i 
is  given  in  terms  of  the  Reynolds  number  R^l  and 


Prandtl  number  Pr  of  the  flow.  In  free  convection, 
where  there  is  no  external  velocity  scale,  Nui^  is 
given  in  terms  of  the  Rayleigh  number  R^l,  which 
is  the  product  of  the  Gr ashof  and  Prandtl  numbers, 
GrL  and  Pr.  Relationships  between  Nwl  and  these 
nondimensional  parameters  have  been  determined 
empirically  for  different  cylinder  roughnesses  and 
orientations. 

To  decide  which  empirical  formula  to  use  in 
forced  convection,  we  need  to  know  the  typical 
range  of  Reynolds  numbers  for  freezing  rain 
accreting  to  a  wire.  This  geometry  is  idealized  as 
flow  around  an  infinite  cylinder  with  the  wire-ice 
diameter  (or  icicle  diameter)  as  the  length  scale. 
The  smallest  diameter  of  interest  is  0,5  cm  for  a 
new  icicle,  and  the  diameter  of  a  moderately  large 
wire-ice  accretion  is  about  12  cm.  Typical  wind 
speeds  during  freezing  rain  range  from  0.5  to  10 
m/s,  and  the  kinematic  viscosity  of  air  is  Va  =  0.132 
cm^/s  at  0°C.  Using  these  values,  the  Reynolds 
number  range  of  interest  is  190  <  R^d  <  9.1  x  10^, 
where  D  is  the  diameter  of  either  the  icicle  or  the 
wire-ice  accretion. 

Achenbach  (1977)  determined  the  total  Nusselt 
number  for  rough  cylinders  in  subcritical  flow,  for 
R^D  >  10^  with  air  {Pr  =  0.72)  as  the  fluid,  as 
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Nun  =  0.18  Rev  (6) 

In  these  wind-tunnel  experiments,  the  cylinder 
surface  was  maiatained  at  a  constant  temperature. 
The  cylinders  were  roughened  to  three  different 
roughness  heights  by  knurling  the  surface  into 
regular  pyramidal  shapes.  The  data  showed  that, 
in  subcritical  flow,  Nud  is  independent  of  the 
roughness  height,  but  the  upper  value  of  the 
Reynolds  number  for  subcritical  flow  varies  be¬ 
tween  2  X 10^  and  10^  decreasing  with  increasing 
roughness.  The  higher  Reynolds  numbers  for  ice 
accretions  in  freezing  rain  are  below  the  upper 
value  of  this  subcritical  flow  range,  but  the  lower 
values  of  Reu  extend  well  below  10^  the  lower 
limit  in  Achenbach's  (1977)  experiments.  To  ob¬ 
tain  appropriate  Nusselt  numbers  for  these  lower 
Reynolds  number  conditions,  I  used  results  from 
Morgan  (1973).  He  determined  Nud  for  bare, 
stranded  conductors  for  Reo  >  10^,  maintaining  a 
constant  heat  flux  rather  than  a  constant  surface 
temperature.  He  compared  these  stranded  con¬ 
ductor  results  to  his  Nusselt  number  formulation 


for  smooth  cylinders  in  Morgan  (1975).  Morgan's 
and  Achenbach's  results  and  a  formula  deter¬ 
mined  for  smooth  cylinders  by  Zhukauskas  and 
recommended  in  Incr opera  and  DeWitt  (1985)  are 
compared  in  Figure  4.  This  figure  shows  that 
Zhukauskas's  smooth  cylinder  formula  predicts 
Nusselt  numbers  that  are  lower  than  Morgan's 
data  at  low  Reynolds  numbers  and  lower  than 
Achenbach's  at  high  Reynolds  numbers.  The  scat¬ 
ter  of  the  wmd-tunnel  data  in  this  figure  illustrates 
the  imcertainty  in  determining  the  Nusselt  num¬ 
ber  for  real  freezing-rain  ice  accretions  with  ir¬ 
regular  roughness  that  is  different  from  that  in 
both  Achenbach's  and  Morgan's  controlled  ex¬ 
periments.  I  incorporated  Morgan's  and 
Achenbach's  rough  cylinder  results  in  the  heat- 
balance  model  as  follows: 

Nuv  =  0.583  Rev  35  <  Rev  <1600 
Nuv  =  0.18  Rev^'^^  1600  <  Keo  <>  1  x  10^  (7) 

Nuv  =  0.00257  Rev>lx  W  , 

extrapolating  Achenbach's  formula  to  lower 


Figure  4.  Nusselt  numbers  for  forced  convection  over  rough  and  smooth  cylinders  as 
a  function  of  Reynolds  number.  Rough  cylinder  data  from  Achenbach  (1977) 
are  for  Reynolds  numbers  larger  than  10^.  Morgan's  (1973)  rough  cylinder  data 
are  for  Reynolds  numbers  less  than  5x10^.  Zhukauskas's  curve  (Incropera  and  DeWitt 
1985)  is  for  smooth  cylinders  for  the  entire  range  of  Reynolds  numbers  shown. 
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Reynolds  numbers  to  intersect  Morgan's  at  Rcd  = 
1600.  In  Morgan's  formula  for  Nud,  Va  and  pa  are 
supposed  to  be  determined  at  the  average  of  the 
wire-ice  surface  temperature  (0°C)  and  the  air 
temperature  T.  However,  these  temperatures  typi¬ 
cally  differ  by  only  a  few  degrees  in  freezing  rain, 
so  there  is  very  little  error  in  using  the  air  tem¬ 
perature  in  calculating  these  parameters. 

Both  Achenbach's  and  Morgan's  experiments 
were  done  in  air,  and  neither  specified  the  depen¬ 
dence  of  Nud  on  the  Prandtl  number.  For  air  at  the 
temperatures  and  pressures  appropriate  for  freez¬ 
ing  rain,  Pr  is  essentially  constant.  Following 
Zhukauskas  (Incropera  and  DeWitt  1985),  I  as¬ 
sume  that  Nwd  is  proportional  to  Thus  eq  7 

can  be  written 


Ap^v  —  PwO  l^^PwT 


=  10" 


Rhej 


273.15  7  +  273.15 


R 


(11) 


where  mw  is  the  molecular  weight  of  water  vapor, 
R  is  the  gas  constant,  Rh  is  the  relative  humidity, 
and  ej  and  Bq  are  the  saturation  vapor  pressures  at 
temperatures  7  and  0°C,  respectively.  This  formu¬ 
lation  for  Qe  is  compared  with  that  in  other  ice 
accretion  models  in  Appendix  A. 

The  empirical  formula  used  in  the  heat-balance 
model  for  the  saturation  vapor  pressure  over  wa¬ 
ter, 

ej  =  6.1121(1.0007  + 

3.46  X  10-6pa)el7.502T/(240.97  +  T)  ^  (12) 


Nwd  =  (0.583/Pr'’'^^)  Pr®'^’' 

35  <  Rex)  <  1600 
Nwd  =  (0.18/Pr°-^^) 

1600  <ReD<  1x10^  (8) 

Nwd  =  (0.00257/ Pr°-2^)  Pr°-^^  Rcd”'^^ 
ReD>lxl0^ 

This  formulation  is  used  with  the  heat  and  mass 
transfer  analogy  in  determining  the  evaporative 
cooling  flux. 

Evaporative  cooling 

Under  typical  freezing-rain  conditions,  the  im¬ 
pinging  water  does  not  freeze  immediately.  There 
is  a  vapor  density  gradient  away  from  the  wire 
because  the  wire-ice  surface  is  warmer  than  the  air 
and  the  saturation  vapor  density  increases  with 
temperature.  This  gradient  is  enhanced  if  the  rela¬ 
tive  humidity  of  the  air  is  less  than  100%.  Thus,  the 
water  film  on  the  surface  of  the  wire-ice  accretion 
is  evaporating  as  it  freezes  and  cooling  the  remain¬ 
ing  water.  The  evaporative  cooling  flux  is 

Qe  =  -10^  ^  =  -10^514  VAPw  '  (9) 

where  Le  is  the  latent  heat  of  evaporation,  me  is  the 
mass  of  evaporated  water  per  unit  length,  is  the 
wire-ice  diameter,  dt  is  the  time  interval,  is  the 
vapor  transfer  coefficient,  and  is  the  vapor 
density  difference  between  the  ice  accretion  sur¬ 
face  and  the  air.  The  vapor  density  pwT  is  related  to 
the  saturation  vapor  pressure  ej  by  the  ideal  gas 
law 

=  10^R(T  +  273.15) ,  (10) 

PwT 

to  give 


is  from  Buck  (1981).  The  saturation  vapor  pressure 
over  water,  rather  than  over  ice,  is  used,  both  in  the 
air  and  at  the  accretion  surface,  because  both  the 
precipitation  and  the  accretion  surface  are  unfro¬ 
zen. 

The  vapor  transfer  coefficient  hm  is  determined 
using  the  mass  and  heat  transfer  analogy  (Incropera 
and  DeWitt  1985).  The  nondimensional  number 
for  mass  transfer  equivalent  to  Nwl  for  heat  trans¬ 
fer  is  the  Sherwood  number  S/zl  =  ^mh/Km  where 
Km  is  the  vapor  diffusivity.  The  Schmidt  number 
Sc  =  Va/Km  is  equivalent  to  Pr.  Like  Pr,  Sc  is 
virtually  constant  for  air  at  the  temperatures  and 
pressures  we  are  concerned  with.  The  heat  and 
mass  transfer  analogy  states  that  for  any  function 
Nwl(P^/  the  Sherwood  number  is  the  same 
function  of  KeL  and  Sc,  Thus,  for  mass  transfer  the 
equivalent  to  eq  8  is 

ShD  =  (0.583/Pr°-^^)  Rcd”'^^^ 

35  <  Rcd  <  1600 

Shv  =  (0.18/Pr°-^^)  ReD°-^^ 

1600  <ReD<  1x10^  (13) 

S/iD  =  (0.00257/Pr“-^^)  Rcd®'^® 

Reo  >  1 X 10^ 

The  vapor  transfer  coefficient  /zm  =  (S/iD^m)/!^  is 
then  obtained  from  eq  13  for  the  appropriate 
Reynolds  number  range.  Significantly  different 
values  for  vapor  diffusivity  Km  are  specified  in 
different  references.  Incropera  and  DeWitt  (1985) 
give  Km  =  0.26  cm^/s  at  7  =  25°C  and  Pa  =  1000 
mbar  in  their  Table  A8.  They  also  say  Km  is  propor¬ 
tional  to  (7  +  273.15)^/^,  which  results  in  Km  =  0.23 
cm^/s  at  0°C.  Batchelor  (1970,  Appendix  1)  speci¬ 
fies  Km  =  0.25  cm^/s  at  15°C.  Pruppacher  and  Klett 
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(1980,  p.  413)  give 


which  results  in  Km  =  0.21  cm  / s  at  Pa  =  1000  mbar 
and  T  =  0°C.  They  also  say  that  conventionally 
used  values  may  be  too  high,  and  as  they  are  a 
better  source  for  this  information  than  either 
Incropera  and  DeWitt  (1985)  or  Batchelor  (1970), 
their  formula  for  Km  is  used  in  the  heat-balance 
model. 

The  mass  of  water  that  evaporates  in  the  pro¬ 
cess  of  freezing  the  impinging  rain  is  mg  = 
-10~^QeDcdt/Le  (eq  9).  In  the  heat-balance  model, 
the  evaporated  water  is  subtracted  from  the  run¬ 
off  water  available  for  icicle  formation  if  not 
all  the  impinging  precipitation  freezes  directly  to 
the  wire. 


Convective  and  evaporative 
cooling  in  no  wind 

When  there  is  no  wind  there  cannot  be  forced 
convection,  so  convective  heat  transfer  occurs  by 
free  convection  only.  Free  convection  occurs  when 
the  temperature  difference  between  the  wire-ice 
accretion  and  the  air  causes  the  air  column  to  be 
unstable.  The  wire-ice  accretion  is  typically  warmer 
than  the  ambient  air,  so  the  warmer  air  over  the 
wire  rises  and  is  replaced  continuously  by  colder 
air.  Heat  transfer  in  free  convection  is  governed  by 
the  Grashof  number, 

^  (15) 

(273.15 +  r)v^ 


and  the  Rayleigh  number. 


P«L  =  GrLPr  = 


g(P-T)L\c^ 
(273.15 +  T)v,fca  ' 


(16) 


where  g  is  the  acceleration  of  gravity,  Cp  is  the 
specific  heat  of  air,  and  Va  is  the  kinematic  viscosity 
of  air.  It  can  be  assumed  that  forced  and  free 
convection  do  not  occur  simultaneously  when 
Reh/Gr\}  is  much  different  from  1  (Incropera  and 
DeWitt  1985).  For  typical  wire-ice  diameters,  wind 
speeds,  and  air  temperatures  in  freezing  rain,  Gr 
«  so  this  is  a  good  assumption.  I  assume  in 
the  heat-balance  model  that  forced  convection 
controls  the  heat  transfer  if  there  is  any  wind,  and 
free  convection  controls  only  when  there  is  no 
wind  {V  =  0).  The  Nusselt  number  determined  by 
Morgan  (1975)  for  free  convection  aroimd  smooth 
horizontal  cylinders  is 


Nmdc  =  0.85  100  <  Raoc  <  10^  (17) 

Nwdc  =  0.48  Ravc°-^°  10^  <  Kadc  <  10^. 

In  a  separate  paper,  Morgan  (1973)  argues  that 
for  relatively  low  Rayleigh  numbers  (100  <Ra  < 
4  X 10^),  the  roughness  of  the  surface  does  not  have 
an  effect  on  heat  transfer,  so  these  formulas  should 
also  apply  to  rough  ice  accretions.  For  most  rea¬ 
sonable  combinations  of  the  wire-ice  accretion 
diameter  and  weather  parameters  in  freezing  rain, 
the  Rayleigh  number  is  within  this  low  range,  so 
eq  17  are  incorporated  in  the  heat-balance  model 
to  determine  the  heat  transfer  coefficient  for  the 
wire  when  there  is  no  wind. 

The  heat  and  mass  transfer  analogy  can  be  used 
in  free  convection  (Incropera  and  DeWitt  1985)  to 
estimate  the  mass  transfer  coefficient  for  evapora¬ 
tive  cooling  from  eq  17: 

S/!Dc  =  0.85  (Sc/Pr)”'^®® 

100  <  Kaoc  <  10^ 

S/iDc  =  0.48  {Sc/Prf^^  (18) 

10^  <  RflDc  <  10^. 

Longwave  radiative  cooling 

In  computing  the  flux  of  longwave  radiation  for 
the  heat-balance  model,  I  assume  that  the  wire-ice 
accretion  is  radiating  at  0°C  and  the  clouds  and 
surrounding  rain-filled  air  are  radiating  at  the  air 
temperature  T.  The  emissivity  of  everything  is 
taken  to  be  1  [emissivity  of  ice  is  0.98  and  of  water 
is  0.96  (Incropera  and  DeWitt  1985)].  Then  the 
longwave  cooling  flux  is  given  by  the  blackbody 
radiation  law: 

Ql  =  -7ca[273.15^  -  (T+273.15)^] .  (19) 

Shortwave  radiative  heating 

The  shortwave  heat  flux  Qs  =  nSr/l  is  deter¬ 
mined  using  the  measured  or  modeled  incoming 
global  diffuse  radiation  flux  Sr.  This  information  is 
available  hourly  at  some  weather  stations  (NOAA 
1993)  and  is  reported  as  an  average  for  each  hour. 
The  factor  of  7c/2  transfers  the  reported  values  for 
a  flat  horizontal  surface  to  a  horizontal  cylinder. 
The  derivation  of  this  factor  is  presented  in  Ap¬ 
pendix  B.  I  assume  in  the  heat-balance  model  that 
the  albedo  of  the  wire-ice  accretion  is  zero  (the  ice 
is  clear),  that  all  the  incoming  radiation  to  the 
accretion  is  absorbed  in  either  the  ice  or  the  under¬ 
lying  wire,  and  that  there  is  no  radiation  reflected 
back  from  the  ground  to  the  wire.  These  assump¬ 
tions  are  only  approximately  correct,  but  they  are 
somewhat  offsetting  in  their  simplicity:  a)  the  al- 
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bedo  of  the  ice  is  close  to,  but  greater  than  zero,  and 
depends  on  how  much  air  and/or  snow  is  incor¬ 
porated  in  the  accretion,  b)  the  shorter- wavelength 
visible  radiation  may  not  be  absorbed  in  the  ice, 
but  it  is  probably  absorbed  in  the  underlying  wire, 
which  is  assumed  to  be  in  thermal  equilibrium 
with  the  ice  accretion,  and  c)  the  ground  surface 
reflects  some  incoming  radiation,  depending  on 
the  ground  cover  and  the  amount  and  condition  of 
snow  and  ice  on  the  groimd. 

The  shortwave  and  longwave  fluxes  are  oppo¬ 
site  in  sign  and  comparable  in  size,  but  there  is  no 
shortwave  flux  at  night,  when  the  air  temperature 
may  be  relatively  cold  and  the  longwave  flux 
relatively  large.  Both  are  small  compared  with  the 
convective  and  evaporative  fluxes  in  typical  freez¬ 
ing-rain  conditions,  but  large  compared  with  the 
viscous  and  kinetic  heat  fluxes  (Fig.  3). 

Viscous  heating 

The  formula  for  the  viscous  heat  flux  is  from 
Makkonen  (1984).  At  wind  speeds  appropriate  for 
freezing  rain,  viscous  heating  from  air  moving  by 
the  wire-ice  accretion  is  small  compared  with  the 
other  heat  fluxes  (Fig.  3).  This  term  is  not  included 
in  the  heat-balance  model. 

Kinetic  heating 

When  precipitation  drops  collide  with  the  wire- 
ice  accretion,  some  of  their  kinetic  energy  is  con¬ 
verted  to  heat.  At  the  droplet  speeds  associated 
with  freezing  rain,  this  term  is  small  (Fig.  3).  It  is 
not  included  in  the  heat-balance  model. 

Resistive  heating 

If  the  wire  conducts  electricity,  heat  is  gener¬ 
ated  from  the  wire  resistance.  The  resistive  heat 
generated  per  unit  length  of  wire  I^Rc  depends  on 
the  current  I  and  the  wire  resistance  per  meter  Rq, 
which  generally  is  a  function  of  the  current  load 
and  frequency.  The  resistive  heat  flux  per  wire 
diameter  is 

Q,(W/m2)  =  15^.  (20) 

The  resistive  heat  flux  in  a  wire  varies  from 
zero,  when  there  is  no  current,  up  to  values  larger 
than  any  of  the  other  heating  or  cooling  fluxes.  For 
example,  for  a  2.63-cm  (1.036-in.)-diameter  ACSR 
conductor  with  I  =  400  amp  and  Rq  =  8.14  x  10”^ 
ohm/m  (0.131  ohms/mile),  Qr  =  490  W/m^.  This 
is  comparable  to  the  typical  convective  and  evapo¬ 
rative  cooling  fluxes  in  freezing  rain  shown  in 
Figure  3.  Because  of  the  extremely  large  range  in 


the  resistive  heating  flux  and  because  assuming  no 
resistive  heating  is  both  reasonable  and  conserva¬ 
tive,  it  is  not  included  in  the  model. 

Latent  heat  of  fusion 

To  freeze  water  at  0°C,  the  heat  of  fusion  must 
be  removed.  The  fraction/ of  the  impinging  pre¬ 
cipitation  that  freezes  is  the  ratio  of  the  sum  of  all 
the  other  heat  fluxes  to  the  latent  heat  flux  Qf  that 
is  calculated  assuming  that  all  the  incoming  water 
freezes: 

r  _  Qc  Qe  Qi  Qs  Qw  ^21) 

^  Qf 

If/ is  1  or  greater,  then  all  the  incoming  water 
freezes.  Properly,  the  temperature  of  the  accretion 
surface  shouldbe  determined  iteratively,  and  heat- 
flux  terms  for  cooling  the  ice  to  the  calculated 
surface  temperature  and  depositing  hoarfrost  or 
sublimating  accreted  ice  should  be  included.  These 
are  small  effects  in  freezing-rain  conditions,  how¬ 
ever,  and  are  not  incorporated  in  the  heat-balance 
model. 

Other  effects 

Two  other,  possibly  significant,  factors  in  de¬ 
termining  the  heat  balance  at  the  ice  accretion 
surface  are  a)  the  angle  of  the  wire  to  the  wind 
direction,  and  b)  freezing  of  the  impinging  pre¬ 
cipitation  by  the  cold  stored  in  the  wire. 

Wind  direction  is  archived  with  the  weather 
data,  so  for  a  particular  structure  the  yaw  of  the 
wire  to  the  wind  direction  during  freezing  rain 
could  be  determined.  The  wire  orientation  could 
then  be  incorporated  in  the  calculation  of  the  pre¬ 
cipitation  flux  zv  across  the  wire  by  multiplying  the 
righthand  side  of  eq  3  by  the  sine  of  the  angle 
between  the  wire  and  the  direction  of  the  precipi¬ 
tation  flux.  The  Nusselt  number  also  depends  on 
the  yaw  of  the  wire.  Morgan  (1973, 1975)  discusses 
the  variation  of  the  Nusselt  number  with  yaw  in 
forced  convection  for  both  smooth  cylinders  and 
bare,  stranded  wires. 

If  the  air  temperature  is  cold  before  freezing 
rain  begins  and  the  wire  is  in  equilibrium  with  the 
air,  impinging  precipitation  may  freeze  by  con¬ 
ductive  heat  loss  into  the  wire.  The  mass  of  water 
that  can  be  frozen  by  this  mechanism  is 

M  = - (22) 

4Lf 

where  pc  is  the  density  and  Cpc  is  the  heat  capacity 
of  the  wire.  This  effect  is  small  except  perhaps  for 
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large  structural  elements.  For  example,  for  a  3-cm- 
diameter  aluminum  wire,  with  Cpc  =  0.88  J/  g°  C,  Pc 
=  2.77  g/cm^,  and  air  temperature  T  =  -S'^C,  the 
mass  of  water  that  can  be  frozen  by  heat  conduc¬ 
tion  into  the  wire  is  26  g/m,  which  is  equivalent  to 
an  ice  thickness  of  0.03  cm.  The  calculation  of  this 
additional  mass  of  ice  is  included  as  an  option  in 
the  model. 


5.  ICICLES 

Precipitation  that  does  not  freeze  to  the  wire 
may  freeze  as  icicles  in  the  process  of  dripping  off 
the  wire.  I  assume  that  there  is  no  precipitation 
impinging  directly  on  the  icicles.  At  low  wind 
speeds,  icicles  form  directly  under  the  wire  and  are 
shielded  by  the  wire  from  the  falling  rain.  At  high 
wind  speeds,  the  runoff  water  is  blown  to  the  lee 
side  of  the  wire  so  icicles  form  behind  the  wire, 
shielded  from  the  wind-driven  rain.  In  the  heat- 
balance  model,  a  single  icicle  per  meter  of  wire 
represents  the  actual  icicles  on  that  meter  of  wire 
stacked  end  to  end.  If  there  is  no  existing  icicle,  the 
mass  of  water  that  is  not  frozen  or  evaporated 
while  on  the  wire  forms  a  protoicicle,  which  is  a 
stack  of  partially  frozen  hemispherical  droplets 
with  an  assumed  diameter  dd  =  0.5  cm  (Maeno  et 
al.  1994).  The  protoicicle  freezes  from  the  outside 
in,  and  when  the  surface  water  freezes  the  unfro¬ 
zen  interior  water  is  retained  and,  as  more  heat  is 
removed,  eventually  freezes.  Assuming  the  sur¬ 
face  layer  of  ice  on  the  hemispherical  droplets  is  0.1 
cm  thick  (Fig.  2,  Maeno  et  al.  1994),  the  ratio  of  the 
volume  of  ice  to  the  volume  of  the  hemisphere  is 
0.64.  In  the  heat-balance  model,  all  the  water  in  the 
protoicicle  is  incorporated  in  the  icicle  if  the  freez¬ 
ing  fraction /b,  determined  from  the  heat  balance 
for  the  protoicicle,  is  greater  than  0.64. 

The  convective  and  evaporative  heat  fluxes  for 
this  protoicicle  are  determined  using  Nudd  for 
spheres.  For  forced  convection  (Incropera  and 
DeWittl985,  eq  7.55), 


tained  from  eq  23  and  24  using  the  heat  and  mass 
transfer  analogy  and  replacing  Pr  with  Sc  both 
explicitly  and  in  Rudd  (Incropera  and  DeWitt  1985). 

Of  the  other  heat-flux  terms,  Qw  =  0  for  icicles 
because  the  runoff  water  is  already  at  0°C,  I  as¬ 
sume  that  there  is  no  solar  radiation  absorbed  in 
the  icicles,  Qs  =  0,  because  a)  the  icicle  orientation 
is  close  to  vertical  and  thus  roughly  parallel  to  the 
incoming  diffuse  solar  radiation  flux,  b)  much  of 
the  longer-wavelength  visible  radiation  is  absorbed 
in  the  cloud  cover,  and  c)  very  little  of  whatever 
incident  shorter-wave  radiation  that  gets  through 
the  clouds  will  be  absorbed  in  the  small-diameter 
clear  icicle. 

The  f r action/b  of  the  available  water  that  freezes 
is  determined  as  in  eq  21.  If/b  is  greater  than  0.64 
then  all  the  available  protoicicle  water  is  retained, 
and  ifyi,  is  less  than  0.64  the  fraction  of  the  protoicicle 
mass  that  is  incorporated  in  the  icicle  is  jb/0.64. 

After  an  icicle  forms,  runoff  water  may  freeze  to 
its  sides.  If  there  is  wind,  the  Nusselt  and  Sherwood 
numbers  for  forced  convection  are  determined 
using  eq  7  and  13,  based  on  the  diameter  of  the 
icicle  Di.  If  there  is  no  wind,  heat  transfer  is  by  free 
convection.  Heat  transfer  by  free  convection  from 
the  icicle  (vertical  cylinder)  is  different  from  that 
for  the  wire  (horizontal  cylinder).  The  relatively 
warm  air  at  the  surface  of  the  icicle  flows  upward, 
drawing  cooler  air  after  it.  The  boundary  layer 
thickness  increases  from  the  bottom  to  the  top  of 
the  icicle,  decreasing  the  heat  transfer.  The  Nusselt 
number  for  free  convection  for  this  geometry  is 
based  on  the  length  of  the  icicles  rather  than  on 
their  diameter.  The  appropriate  icicle  length  to  use 
is  that  for  a  single  icicle,  rather  than  this  stack  of 
icicles.  Following  Makkonen  and  Fujii  (1993),  I 
assume  45  identical  icicles  per  meter  of  wire  to  give 
an  average  icicle  length  Lbar  =  Pi/ 45.  Because  the 
icicles  are  closely  spaced  on  the  wire,  thus  inhibit¬ 
ing  cooling,  I  use  an  empirical  formula  for  heat 
transfer  from  a  vertical  plate,  rather  than  from  a 
cylinder  to  determine  the  Nusselt  number 
(Incropera  and  DeWitt  1985,  eq  9.27): 


Nudd  =  2  +  Pr°  \0A  Redd^^^ 

+  0.06  Redd^/%  (23) 


and  for  free  convection  (Incropera  and  DeWitt 
1985,  eq  9.35), 


0.589RflH 


1+{0A69/Prf^^^f 


The  Sherwood  number  for  vapor  transfer  is  ob¬ 


The  Sherwood  number  for  evaporative  cooling  is 
obtained  from  eq  25  by  replacing  Pr  with  Sc,  both 
explicitly  and  in  R^Lbar-  The  convective  and  evapo¬ 
rative  heat  fluxes  are  quite  insensitive  to  the  as- 
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S/D  =  7t 


S/D  =  2.5 


S/D  =  4.0 


0.35  cm 


0.28  cm 


0.44  cm 


Figure  5.  Equivalent  uniform  radial  ice  accretions  on  cylinders  with  circular,  rectangular,  and  angle 
cross  sections  for  a  rainfall  ofl  cm.  The  width  of  the  rectangle  is  1/4  of  its  length.  The  angle  has  equal- 
length  legs  with  arbitrary  thickness.  The  perimeter-to-"  diameter"  ratio  S/D  ranges  from  3.14  for  the 
circle  to  2.5  for  the  rectangle  and  4  for  the  angle. 


sumed  average  icicle  length  as  both  h  and  hm  are 
proportional  to  Lbar“^/^-  As  for  the  protoicicle,  Qw 
=  0  and  Qs  =  0  and  the  freezing  fraction  fs  is 
determined  using  eq  21.  The  water  that  does  not 
evaporate  or  freeze  to  the  sides  of  the  icicle  forms 
a  protoicicle  at  its  tip,  and  the  heat  balance  for  the 
protoicicle  is  determined  as  before.  Finally,  the 
water  that  does  not  freeze  at  the  tip  is  assumed  to 
drip  off.  The  new  total  icicle  mass  Mj  is  distributed 
evenly  over  a  cylinder  with  length  Li  equal  to  the 
original  icicle  length  plus  the  length  of  the 
protoicicle. 


6.  SIMPLE  FLUX  MODEL 

The  simplest  case  of  ice  accreting  in  freezing 
rain  occurs  if  all  the  precipitation  impinging  on  the 
structure  freezes  in  a  imiform  radial  accretion. 
When  that  happens,  the  amount  of  ice  on  the 
structure  is  directly  related  to  the  amount  of  ram 
that  falls.  For  example,  if  1  cm  of  freezing  rain  falls 
on  a  2“Cm-wide  flat  plate  and  freezes  in  place,  then 
a  layer  of  ice  l(po/pi)  =  1.1  cm  thick  results.  On  a  4- 
cm-wide  plate  nearby,  the  same  thickness  of  ice 
would  form  with  twice  the  mass.  A  2-cm-diameter 
horizontal  circular  cylinder  intercepts  the  same 
depth  of  rain  as  the  2-cm-wide  flat  plate.  If  that 
water  depth  is  then  spread  imiformly  around  the 
cylinder's  circumference  and  frozen,  forming  a 
imiform  radial  accretion,  then  the  layer  of  ice  is 
l(Po/pi)/7r  =  0.35  cm  thick.  The  factor  of  n  is  the 
ratio  between  the  circular  cylinder's  circumfer¬ 
ence,  over  which  the  ice  freezes,  and  the  diameter. 


which  intercepts  the  precipitation.  This  calcula¬ 
tion  can  be  generalized  to  a  cylinder  with  any  cross 
section  to  determine  the  uniform  ice  thickness  on 
its  perimeter: 


Rgq  — 


SPi  ;=1 


(26) 


where  the  summation  gives  the  total  precipitation 
amount  in  N  hours  and  S  and  D  are  the  perimeter 
and  horizontal  dimension,  respectively,  of  the  cross 
section.  Equivalent  uniform  ice  thicknesses  are 
shown  in  Figure  5  for  three  different  cross  sections 
in  a  1-cm  rainfall. 

For  many  shapes,  but  not  circles,  the  D/S  ratio 
decreases  as  ice  accretes  imiformly  on  the  perim¬ 
eter.  For  those,  Req  could  be  determined  incremen¬ 
tally.  When  this  is  done  for  the  rectangle  in  Figure 
5,  using  0.2-cm  increments  of  rain,  the  final  uni¬ 
form  radial  ice  thickness  is  0.42  cm  rather  than  0.44 
cm.  This  difference  is  well  within  the  uncertainties 
in  either  this  simple  model  or  the  detailed  heat- 
balance  model. 

This  simple  precipitation  model  shows  that,  at 
least  to  first  order,  the  thicknesses  of  ice  on  compo¬ 
nents  with  the  same  shape  cross-section  but  differ¬ 
ent  diameters  are  the  same.  On  components  with 
the  same  diameter  but  different  cross  sections,  the 
ice  thickness  ratio  is  equal  to  the  inverse  of  the 
perimeter  ratio.  Thus,  uniform  radial  ice  thick¬ 
nesses  determined  by  this  model,  or  any  other 
model,  for  a  circular  cylinder  of  any  diameter  can 
be  extrapolated  both  to  cylinders  with  other  cross 
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sections  and  cylinders  with  other  diameters. 

A  slightly  more  complex  model  takes  the  in¬ 
creased  flux  of  water  due  to  the  horizontal  velocity 
of  the  rain  drops  in  the  wind  into  accoxmt.  Using 
the  flux  of  precipitation  from  eq  3,  the  uniform 
radial  ice  thickness  is 


R 


eq  ” 


(27) 


This  simple  flux  model  with  S/D  =  n  is  incorpo¬ 
rated  in  ZRAIN  for  comparison  with  the  heat- 
balance  model. 


7.  COLLISION  EFFICIENCY 

In  in-cloud  icing,  the  collision  efficiency  of  the 
cloud  droplets  with  a  structure  is  an  important 
factor  in  determining  the  rate  of  ice  accretion. 
Cloud  droplets  are  very  small,  typically  5  to  50  pm 
in  diameter,  and  have  correspondingly  small  ter¬ 
minal  velocities  (0.007  m/s  for  15  pm  droplets. 
Best  1950).  Thus,  in  windless  conditions  there  is 
essentially  no  droplet  flux.  When  there  is  wind,  the 
wind  drag  on  the  droplets  carries  them  along  and 
the  droplets  follow  the  wind  streamlines  aroimd 
any  obstacle  in  their  path.  Only  the  inertia  of  the 
droplets,  which  tends  to  keep  them  moving  in  a 
straight-line  path,  makes  them  diverge  from  the 
wind  streamlines.  The  collision  efficiency  of  the 
cloud  droplets  with  a  structure  represents  a  bal¬ 
ance  between  their  drag,  the  tendency  to  follow 
the  streamlines,  and  their  inertia,  the  tendency  to 
continue  in  a  straight  line.  Qualitatively,  the  re¬ 
sults  of  this  balance  are  that  a)  smaller  droplets, 
which  have  less  inertia,  have  lower  collision  effi¬ 
ciencies,  b)  the  collision  efficiency  of  droplets  is 
smaller  for  large  obstacles  than  for  small  ones, 
because  streamlines  diverge  relatively  farther  in 
front  of  large  obstacles  than  small  ones,  and  c)  in 
high  winds  droplets  have  more  momentum  so 
their  collision  efficiency  is  higher.  These  observa¬ 
tions  are  quantified  in  Langmuir  and  Blodgett 
(1946).  They  numerically  solved  the  trajectory  equa¬ 
tion,  ignoring  gravitational  effects,  to  determine 
droplet  collision  efficiency  as  a  function  of  two 
nondimensional  numbers  that  depend  on  wind 
speed,  wire  diameter,  droplet  diameter,  and  the 
density  and  viscosity  of  air. 

The  Langmuir  and  Blodgett  results  cannot  be 
applied  to  precipitation  droplets.  The  terminal 


velocity  of  drizzle  and  rain  drops,  typically  100s  of 
microns  to  1  or  2  mm  in  diameter,  is  between  about 
2  and  7  m/s  (Wang  and  Pruppacher  1977).  Thus, 
even  in  no  wind,  there  is  a  downward  flux  of  water 
that  is  a  balance  between  gravity  and  the  drag  of 
the  still  air  on  the  drops.  The  drops  fall  through  the 
air  and  collide  either  with  the  ground  or  with 
structures  in  their  path.  When  there  is  wind,  the 
drops  gain  horizontal  velocity  in  addition  to  their 
vertical  velocity.  The  effect  of  wind  on  drizzle-  and 
rain-sized  drops  could  be  determined  by  numeri¬ 
cally  integrating  the  trajectory  equation  (Lozo wski 
and  Oleskiw  1983),  including  the  gravitational 
term: 


dt  MPq 


v-V  (5-Pj  +  f  / 


(28) 


where  v  is  the  drop  velocity,  V  is  the  wind  veloc¬ 
ity,  Cd  is  the  drag  coefficient  of  the  drops,  d  is  the 
drop  diameter,  and  g  is  the  acceleration  of  grav¬ 
ity.  In  the  absence  of  a  solution  to  this  equation  for 
rain-  and  drizzle-sized  drops,  it  is  reasonable  to 
assume  that  the  drops  fall  at  their  terminal  veloc¬ 
ity,  move  horizontally  at  the  wind  speed,  and 
collide  with  all  obstacles  in  their  path. 

In  some  freezing-rain  ice  accretion  models,  the 
collision  efficiency  of  the  rain  and  drizzle  drops  is 
calculated  incorrectly.  The  resultant  velocity  of  the 
wind  speed  and  the  drops'  terminal  velocity  is 
used  in  the  Langmuir  and  Blodgett  (1946)  formu¬ 
lation  to  determine  the  collision  efficiency.  This 
ignores  the  different  physics  of  falling  and  wind- 
carried  droplets  and  leads  to  obviously  ridiculous 
results:  0.5-mm-diameter  drizzle  drops  falling  at  2 
m/s  in  no  wind  are  calculated  to  have  a  collision 
efficiency  of  0.56  with  a  1-m-diameter  cylinder 
and  0.05  with  a  10-m  cylinder.  In  that  world,  if  you 
carried  a  big  enough  umbrella,  not  only  would 
you  stay  dry,  so  would  the  umbrella. 


8.  WIND  LOAD 

It  is  often  important  to  know  the  wind  load  on 
a  structure  both  during  a  freezing-rain  storm,  and 
for  as  long  after  the  storm  as  ice  remains  on  the 
structure.  The  projected  area  of  the  structure  is 
larger  because  of  the  ice  accretion,  so  at  a  given 
wind  speed  the  wind  load  is  greater  than  it  would 
be  on  the  bare  structure.  In  ZRAIN,  the  wind  load 
per  meter  of  wire  is  determined  every  hour  using 
the  wire-ice  diameter  and  the  icicle  diameter  and 
length,  assuming  a  drag  coefficient  Cd  =  1  for  both: 
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^load=0.5p,y2[D,  +  Di^].  (29) 

The  factor  Li  / 100  converts  the  wind  load  per  length 
of  icicle  to  wind  load  per  meter  of  wire.  The  wind 
load  results  are  useful  for  identifying  the  combina¬ 
tion  of  wind  and  ice  in  each  event  that  causes  the 
highest  horizontal  load.  This  combination  is  inde¬ 
pendent  of  drag  coefficient  as  long  as  it  can  be 
assumed  to  be  the  same  for  both  the  wire-ice 
accretion  and  the  icicle. 


9.  SAMPLE  RESULTS 

ZRAIN  was  lun  for  the  freezing-rain  events 
that  occurred  at  the  Des  Moines,  Iowa,  airport 
between  1948  and  1993.  The  equivalent  radial  ice 
thicknesses  on  a  2.54-cm-diameter  wire  10  m  above 
the  ground,  hindcast  using  the  heat-balance  model, 
are  shown  in  Figure  6  for  the  316  freezing-rain 
events  in  those  45  years.  For  many  events  either  no 
ice  or  very  little  ice  accreted,  even  though  freezing 
rain  or  drizzle  was  observed,  typically  because 
there  was  little  or  no  measured  precipitation.  The 
largest  accretion  predicted  by  the  model  occurred 
in  January  1982.  That  event  consisted  of  two  freez¬ 
ing-rain  storms  separated  by  19  days  of  cold 
weather  with  temperatures  as  low  as  -30°C.  In  the 
first  12  hours  of  freezing  rain,  0.75  cm  of  ice  accreted, 
and  in  the  second  19-hr  period  of  freezing  rain 
there  was  an  additional  0.75-cm  uniform  radial 
accretion. 

It  is  interesting  to  compare  the  modeled  ice 
thicknesses  with  the  Storm  Data  descriptions  of  the 
severe  ice  storms  in  1990  and  1991.  For  the  7-8 
March  1990  storm  at  the  Des  Moines  airport,  the 


heat-balance  model  hindcast  a  0.5-cm-thick  uni¬ 
form  radial  accretion  with  26-cm-long  icicles  (as¬ 
suming  45  icicles /m),  which  is  equivalent  to  a  1 .3- 
cm  imif orm  radial  ice  accretion.  According  to  S  torm 
Data  this  ice  storm  affected  most  of  Iowa  and 
resulted  in  1.3  to  7.6  cm  (0.5  to  3  in.)  of  ice  on 
exposed  surfaces  and  7.6-  to  12.7-cm  (3-  to  5-in.)- 
thick  ice  on  some  power  lines.  Central  Iowa  from 
near  Des  Moines  and  up  to  the  north  and  west, 
particularly  Carroll  and  Crawford  counties  (about 
80  miles  from  Des  Moines'  Polk  County),  was 
hardest  hit. 

According  to  the  model,  the  31  October  1991  ice 
storm  was  less  severe  than  the  1990  storm  in  the 
Des  Moines  area.  The  hindcast  uniform  radial  ice 
accretion  is  only  0.2  cm  thick  with  3-cm-long  icicles. 
Storm  Data  describes  this  ice  storm  as  affecting 
about  a  third  of  Iowa,  including  the  forecast  zone 
where  Des  Moines  is  located.  From  2.5  to  5.1  cm  (1 
to  2  in.)  of  ice  accumulation  was  reported  "in  many 
areas,"  and  much  of  the  damage  was  attributed  to 
the  high  winds  that  accompanied  the  storm.  How¬ 
ever,  Des  Moines'  Polk  County  was  not  included 
in  the  Presidential  Disaster  Declaration  for  43  Iowa 
counties  following  that  storm  (FEMA  1991),  and 
the  measured  wind  speeds  at  Des  Moines  airport 
were  lower  during  this  storm  than  in  the  1990 
storm. 

The  simple  flux  model  results  are  compared 
with  the  heat-balance  model  results  in  Figure  7. 
The  models  agree  well  for  the  vast  majority  of  the 
freezing-rain  events  at  Des  Moines .  In  the  cases  for 
which  the  results  from  the  simple  model  are  sig¬ 
nificantly  higher  than  those  from  the  heat-balance 
model,  the  air  temperature  was  at  or  above  0°C 
with  relative  humidity  near  100%,  and  moderate 


Figure  6.  Uniform  radial  ice  thicknesses  hindcast  by  the  heat-balance  model  for  freezing-rain  events  at  the  Des  Moines 
airport  from  1948  to  1993. 
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or  heavy  precipitation.  A  time  series  of  one  of  these 
events,  the  7-8  March  1990  ice  storm,  is  shown  in 
Figure  8.  This  storm  was  characterized  by  high 
precipitation  rates,  up  to  1  cm/hr,  increasing  tem¬ 
perature  and  relative  humidity,  and  high  diffuse 
solar  radiation  during  the  day.  As  the  air  tempera¬ 
ture  increased  and  freezing  rain  continued  near 
the  end  of  Julian  day  66,  little  additionalice  accreted 
directly  on  the  wire,  but  much  of  the  rimoff  water 
froze  as  icicles.  It  is  during  this  time  period  that  the 
heat-balance  model  results  diverge  from  the  simple 
flux  model,  which  accretes  all  the  impinging  preci¬ 
pitation. 

The  heat-balance  model  was  also  run  for  a  10- 
cm-diameter  wire  to  determine  whether  the  lack 
of  dependence  of  the  ice  thickness  on  wire  diam¬ 
eter  predicted  by  the  simple  model  also  holds 
when  the  heat  fluxes,  which  are  diameter-depen- 
dent  (Fig.  3),  are  incorporated  in  the  ice  accretion 
calculation.  The  equivalent  radial  thicknesses  for 
these  two  diameters  are  compared  in  Figure  9  and 
show  almost  no  difference.  This  is  undoubtedly 
due  to  the  accretion  of  runoff  water  as  icicles, 
which  freeze  independent  of  the  wire  diameter. 

One  small  but  well-documented  freezing-rain 
event  at  CRREL  on  28  February  1995  serves  as  a 
test  for  the  two  models  against  reality.  This  event 
was  characterized  by  cold  temperatures,  very  light 
winds,  and  light  freezing  drizzle  that  froze  quickly 


rather  than  as  icicles.  The  ice  thickness  determined 
from  the  measured  mass  of  ice  on  a  2.65-cm- 
diameter  horizontal  rod  (Itagaki,  in  prep.)  at 
CRREL's  freezing-rain  weather  station  was  0.19 
cm.  Both  models  predict  equivalent  radial  ice  ac¬ 
cretions  0.18-cm  thick.  While  this  good  agreement 
is  encouraging,  further  testing  of  the  models  in 
more  extreme  freezing-rain  conditions  is  desirable. 

10.  DISCUSSION 

This  report  describes  two  models  to  calculate 
freezing-rain  ice  loads  onhorizontal  cylinders  from 
weather  data.  One  model  is  very  simple  and  re¬ 
quires  only  the  information  that  freezing  rain  is 
occurring  along  with  the  concurrent  precipitation 
rate  and  wind  speed.  The  more  detailed  model 
requires  as  additional  information  the  air  tem¬ 
perature,  dew-point  temperature,  and  solar  radia¬ 
tion  flux,  as  well  as  the  diameter  of  the  wire  on 
which  ice  is  accreting.  A  comparison  of  the  two 
models  using  the  historical  weather  data  at  Des 
Moines  indicates  that  the  simple  flux  model  is 
conservative,  as  expected,  but  generally  agrees 
well  with  the  detailed  heat-balance  model. 

As  freezing-rain  events  occur  at  CRREL,  mea¬ 
surements  of  accreted  ice  on  horizontal  rods  will 
be  compared  with  the  models'  predictions.  This 


Figure  7.  Comparison  of  uniform  radial  ice  thicknesses  at  Des 
Moines  from  the  simple  flux  model  and  the  heat-balance  model. 
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Figure  8.  Time  series  of  the  weather  conditions  and  hindcast  ice  accretions 
for  the  ice  storm  of  7~8  March  1990  at  the  Des  Moines  airport.  The 
formation  of  icicles  in  the  heat-balance  model  occurs  in  relatively  warm, 
humid  conditions.  In  those  conditions  the  results  from  the  heat-balance 
model  and  the  simple  flux  model  diverge. 


field  data  can  be  used  to  fine-tune  the  heat-balance 
model.  Parameters  that  may  have  a  significant 
effect  on  the  calculated  ice  load  and  may  need  to  be 
adjusted,  include  a)  the  0.64  threshold  for  incorpo¬ 
rating  all  of  the  runoff  water  in  the  protoicicle,  and 
b)  the  projected  area  of  the  ice  accretion  that  is 
assumed  to  intercept  the  precipitation  flux. 

The  good  agreement  between  the  heat-balance 
model  and  the  simple  flux  model  is  interesting.  It 


indicates  that  the  modeled  ice  thicknesses  can  be 
easily  extrapolated  so  that  consistent  ice  thick¬ 
nesses  can  be  used  in  the  design  of  structural 
elements  with  different  cross  sections.  It  also  shows 
that  weather  data  can  be  used  to  determine  conser¬ 
vative,  but  apparently  realistic,  ice  loads  using 
back-of-the-envelope  calculations.  It  will  be  inter¬ 
esting  to  see  if  this  good  agreement  holds  as  the 
heat-balance  model  is  fine-tuned. 
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Figure  9.  Comparison  of  uniform  radial  ice  accretion  thicknesses 
for  2.54-  and  1 0-cm-diameter  wires  using  the  heat-balance  model 
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APPENDIX  A:  COMPARISON  OF  EVAPORATIVE 
HEAT  FLUX  FORMULATIONS 


The  formulation  of  the  evaporative  heat  flux  term  in  the  heat-balance  model  is 
different  from  that  of  Lozowski  et  aL  (1987)  and  Makkonen  (1984).  The  different  QeS 
are  compared  in  this  appendix.  Assuming  consistent  units  and  considering  the  heat 
flux  across  an  arbitrary  surface,  Qe  is  (eq  4) 


Qe  “ 


1,273.15 


RHc-y 

r  +  273.15j'X 


(Al) 


Using  the  definition  of  the  Sherwood  number  S/zd  and  the  relationship  between  S/zd 
and  Ihe  Nusselt  number  Nwd  specified  by  the  heat  and  mass  transfer  analogy  gives 


NuD  =  /!D/fca  =  CPr°-^^ReD*’  (A2) 

Shu  -  /JmD/Kjn  =  CSc°-3^  Reo’’ 


where  C  and  b  depend  on  Reo*  Equation  A2  can  be  used  to  write  /zm  in  terms  of  Nud: 
nO.37 


(A3) 


Using  this  relationship  and  Sc/Pr  =  fca/(Pa^^pKm)/  eq  Al  can  be  written  finally  as 


Q. 


_  Nzzp 
D 


nO.37 


0.63 

""  R 


f  gp _ ^be-j  ^ 

’U73.I5  T  + 273.15  j 


(A4) 


The  evaporative  heat  transfer  term  in  Lozowski  et  al.  (1987)  is  given  as 


Qe=h 


/p,N0.63 
h\  —  eLg 

^ Cq  -Rhej^ 

UcJ 

c^Pa 

\  P  J 

(A5) 


where  e  =  is  the  ratio  of  the  molecular  weights  of  water  vapor  and  air.  Using 

the  ideal  gas  law, 

P^^P,R(T-h  273.15)^  (A6) 

nJa 

and  the  definition  of  Nud  in  eq  A2,  eq  A5  becomes 


Qe  = 


Nup 


D 


nO.37 


PaCp 


0.63 


m^Lei^o  - 
R(r  + 273.15) 


(AT) 


Thus,  Lozowski  and  Gates's  formulation  (eq  A7)  is  different  from  this  model's 
formulation  (eq  A4)  only  in  the  temperature  used  to  calculate  the  saturation  vapor 
density  at  the  accretion  surface. 

In  Makkonen  (1984)  the  evaporative  cooling  flux  is  given  as 


Qe  = 


htLe(eo-ej)  _ 
CpPa 


(A8) 
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Again  using  the  ideal  gas  law  and  eq  A2,  this  can  be  rewritten  as 


N»d  fca  ^t)  (A9) 

®  D  [paCpJ  R(r  + 273.15)  ' 

Makkonen's  formulation  differs  from  this  model's  by  the  temperature  used  to 
calculate  the  saturation  vapor  density  at  the  freezing  surface,  by  the  dependence  on 
fca/(paCp)  rather  than  on  [fca/(paCp)]°’^^Km°'^^  and  by  using  the  saturation  vapor 
density  of  the  air  rather  than  the  vapor  density  at  the  ambient  humidity. 

The  effect  of  these  differences  is  to  make  Lozowski  and  Gates's  Qe  4%  larger  and 
Makkonen's  Qe  32%  smaller  than  this  model's  at  typical  freezing-rain  conditions 
T = -3°C,  Rh  =  90%  and  Pa  =  1000  mbar .  The  differences  between  the  models  increase 
somewhat  as  air  temperature  decreases  and  increase  substantially  as  humidity 
decreases. 

This  comparison  addresses  the  differences  in  the  physical  models  of  evaporative 
cooling  and  in  the  assumed  dependence  of  the  Nusselt  number  on  the  Prandtl  and 
Schmidt  numbers.  There  are  further  differences  between  the  models  in  the  empirical 
formulas  used  to  relate  the  Nusselt  number  to  the  Reynolds  number  and  Rayleigh 
number. 
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APPENDIX  B:  SKY  DIFFUSE  SOLAR  RADIATION  FLUX  TO  A 
HORIZONTAL  CYLINDER 

Sky  diffuse  radiation  under  an  overcast  that  would  be  typical  during  freezing  rain 
is  approximately  isotropic  over  the  sky  dome.  The  flux  of  diffuse  solar  radiation  to 
a  surface  inclined  at  an  angle  a  to  the  horizontal  (Fig.  Bl)  is 

Qst  =  Sr(l  +  cosa)/2  (Bl) 

for  the  top  of  the  surface,  and 

Qsb  =  Sr(l  -  cosa)/2  (B2) 

for  the  bottom  of  the  inclined  surface  (Iqbal  1983).  The  circumference  of  a  cylinder 
is  made  up  of  elements  with  area  dC  =  Dda/2  per  imit  length,  inclined  at  angles  to 
the  horizontal  varying  from  -n/2  to  nf2  (Fig.  B2).  The  contributions  of  the  sky  diffuse 
solar  radiation  to  each  of  these  surface  elements  are  summed,  taking  advantage  of 
the  symmetry  of  the  cylinder  about  the  vertical,  to  determine  the  solar  power  per 
imit  length  of  the  cylinder: 


Ps 


yDA^I^p  +  cosa 

_  ItJ  J  I  2 


1 


-cos  a 


(B3) 


Evaluating  the  integrals  results  in  Pg  =  nDSr/2,  so  the  flux  of  sky  diffuse  solar 
radiation  to  the  cylinder  per  diameter  is 


Qs  = 


nSr 


which  was  given  in  eq  4. 


(B4) 


a  =  K/2 


a=:-7i;/2 


Figure  Bl.  Sky  diffuse  solar  radia-  Figure  B2.  Cylinder  surface  elements, 

tion  incident  on  surface  inclined  at 
an  angle  a  to  the  horizontal. 
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